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Q ■ Abstract 

^ ■ The discrete Schrodinger equation with a quasiperiodic dichotomous poten- 

tial specified by the Fibonacci sequence is known to have a singular continuous 
. , eigenvalue spectrum with all states being critically localized. This equation 

l/-^ , can be transformed into a quasiperiodic skew product dynamical system. In 

§: 

critically localized states correspond to fractal attractors which have all Lya- 
punov exponents equal to zero. This provides an alternate means of studying 
the spectrum, as has been done earlier for the Harper equation. We study 



> 

^ , the spectrum of the Fibonacci system and describe the scaling of gap widths 

U '■ 

^ - with potential strength. 



I. INTRODUCTION 

The phenomenon of localization, first discussed in the context of disordered systems by 
Anderson [|1|, has been of great interest in condensed matter physics for several decades. 
The Anderson model is the tight-binding Hamiltonian 
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where c| is the creation operator at the lattice site indexed by i, and Wi and lij are on-site 
and hopping matrix elements. In the original model, 1^ = I for i,j nearest neighbours and 
zero otherwise, and Wi is chosen randomly in the interval [—W/2,W/2]. I and W can be 
combined into a dimensionless ratio which is a measure of the disorder. For ^ > [^]c, a 
critical value, all the states are exponentially localized. Below the critical value, there is 
a mobility edge and states can be extended. For ordered systems, such as the analogous 
system with a periodic potential, states are extended and the energy spectrum is a band 0. 

There are systems that lie between these extremes, for instance crystals with additional 
periodic modulation which is incommensurate with the period of the underlying lattice . 
The lack of strict translational order has an effect similar to disorder, and thus it is of interest 
to know the nature of the states in such systems. In this paper we study the spectrum of a 
1-d chain, the potential at each site being determined by a quasiperiodic sequence deriving 
from the Fibonacci numbers P-[T^. 

The paradigm for studies of quasiperiodic Schrodinger problems has been the Harper 
equation ||13| 

Ipn+l + V'n-l + Vnlpn = Elpn, (2) 

which is the discrete Schrodinger equation for a particle in the quasiperiodic potential Vn = 
2ecos27r(?T,u; + (po), with u being irrational. This is a tight-binding model describing the 
motion of a electron in two dimensional lattice in the presence of a magnetic field, but 
is also known to arise in numerous other contexts. As a result, it has been extensively 



studied both in the physics [[l4|-pl]l as well as in the mathematics p2|-p6| literature (as 
the "almost-Mathieu" equation). Since the results that are known for this system have 
considerable bearing on the present work, we first briefly review the salient features of the 
Harper system. 

Transfer matrices provide a powerful methodology for the calculation of the band spec- 
trum in matrix form, the discrete Schrodinger equation, Eq.(|D can be written as 



T 

-J- r) 



(3) 



where the transfer matrix T„ is given by 

f E~Vn -l\ 

1 



T 

71. 
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(4) 



For 1-d chains, the state vector at the final site is obtained from the starting value by matrix 
multiplication, 



T 



(5) 



with T = Tjq ■ Tjv„i ■ T/v_2 ■ ■ - Ti. The eigenvalues of T„ can be calculated by the Cayley- 
Hamilton theorem as 



A 



E-Vn±\KE-VnY-A 



(6) 



From this equation it is evident that A can be complex if | — |< 2: this corresponds to 
extended states. \i \ E — Vn \> 2, then A is real, corresponding to localized states. There is 
a metal-insulator transition at the so-called mobility edge \ E — Vn |=2. 

In order to determine the actual nature of the spectrum for irrational c<j, a barrage of 
mathematical techniques have been applied to the study of the Harper equation [p!^ ,^8|-|32 



One of the main results has been the discovery of duality by Andre and Aubry IQ. Since 
the Harper potential contains no harmonics, there can be a solution to Eq. @) such that 



ipn = exp{ifc?7,} ^ exp{im(27rci;n + (^q)} 

m=—oo 

where the Fourier coefficients fm are given by 

CXD 

/m = exp{— ifcm} ^ T/^n exp{— m(27rc<jm + 0o)}- 

m=— oo 

Substituting the value of ipn in the above equation yields 



(7) 



(8) 
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2 1 

fm+i + fm-1 + - cos{2num + k)f^ = -Ef^, (9) 

which, written as 

fm+i + fm-i + 2e cos(27ru;m + /c)/„ = Efm (10) 
is similar in form to Eq. (EI), with e = 1/e and E = E/e. The equations are identical when 



1 



with ip = f. Now if fm is locahzed then the sum Y^^=_oo I /m P is finite, and therefore the 
function 

oo 

f {x) = exTp{ikx} ^ /m exp{im27ra;a;} (11) 

m=— oo 

converges. Thus f{x) behaves hke a Bloch function and represents an extended state. For 
X = n + f{x) = exp{ik-^}iljn- Hence if ipn is extended then J2n I V^n P diverges. 
Therefore, if ipn is extended then is locahzed. 

The duality of the model is this correspondence which maps the region e > 1 to the 
region e < 1, with the point e = 1 being "self dual". Thus states for this value of e are 
neither extended nor localized: they are critical. 

This criticality has two manifestations. The spectrum of the Harper equation at e = 
1 is singular continuous: the eigenvalues form a Cantor set. The wave-functions for these 
critical states are power-law localized. As a function of uj, the spectrum of the Harper 
equation, which was first studied in some detail by Hofstadter has a remarkable shape. 
For every rational value, it consists of a finite number of bands, while for irrational values 
of uj it consists of an apparently fractal set of points, giving what is termed the Hofstadter 
butterfly. 

A number of other quasiperiodic potentials |T6|-p^ have since been shown to support 



critical states. One widely studied system is the Kohmoto model, defined at lattice site n 
through the rule |T6|,p 



Vn = a < {nuj} < oj 

= -a Lu < {nuj} < 1, (12) 

(the notation is {y} = y mod 1,) which is quasiperiodic if u is an irrational number. For 
the case of a; = 7 = {y/E — l)/2, the golden mean ratio, one obtains the so-called Fibonacci 
chain. 

This latter potential is the subject of the present paper. Our approach is based on the 
equivalence between the discrete Schrodinger equation and a derived iterative mapping for 
the amplitude ratio of the wave function at neighbouring sites. This equivalence was first 



noted by Bondeson et al. who showed that the Schrodinger equation with a quasiperiodic 
potential could be transformed into a dynamical system with quasiperiodic forcing. Ketoja 
and Satija ||29| extended this analysis to the Harper equation, Eq. (|]), obtaining the entirely 
equivalent Harper map, 

Xk — + 2e cos 27r0fc 
= + t^}, (14) 

using the transformation ifji-i/ipi Xi. This is a skew-product driven iterative mapping of 
the infinite strip (00, 00) (g) [0, 1] to itself. If the frequency uj is an irrational number, the 
driving is quasiperiodic in time. For the Kohmoto model, the analogous corresponding map 
is 

with Vk given by Eq. ([1^). In either case, the quantum problem is meaningful only when 
E is an eigenvalue, but in the above map, E appears only as a parameter. Such driven 
maps, which were first introduced by Grebogi et al. |^, have been studied extensively in the 
context of strange nonchaotic dynamics PB[. In quasiperiodically driven dissipative systems. 



the dynamics can (for appropriate parameter values) be on fractal sets which concurrently 
have non-positive Lyapunov exponents: such motion is on strange nonchaotic attractors 
(SNAs). 



As has now been demonstrated in several studies, localized states correspond to SNAs 
TT| . The transition from extended to localized states as a function of potential strength 
can be viewed as a transition to SNAs. Critical states have a special significance: these are 
SNAs on which all Lyapunov exponents are zero. In this paper we use this fact to study the 
states of the Fibonacci chain for which it is known explicitly that all states are critical, by 
determining the eigenvalue spectrum of the system through the condition that the Lyapunov 
exponent for the specified value of E be zero. This provides a simpler alternative than the 
transfer matrix techniques, and in addition, gives a characterization of the states of the 



system which is not accessible to purely quantum mechanical methods ||31|| . 

In the following section we briefly review the relevant features of SNAs and the methods 
used to study them. Section III contains the main results of this work for the Fibonacci 
chain. Section IV contains a discussion of related potentials deriving from abstract aperiodic 
sequences, followed by a brief summary. 

II. STRANGE NONCHAOTIC ATTRACTORS 

Strange nonchaotic attractors are frequently found in nonlinear systems where the forcing 



is quasiperiodic [|36|. There are similarities between SNAs and both periodic as well as 
chaotic attractors. Like the former, they are characterized by zero or nonpositive Lyapunov 
exponents, and like the latter, they have a fractal structure |^5|]. Owing to fractal structure 
of SNAs, the dynamics is strictly aperiodic. Numerous examples are known now of systems 
with SNAs. These include, in addition to discrete quasiperiodically forced maps, a number 
of continuous dynamical systems such as driven pendulums and oscillators ||36|j37|] . 



In the mapping first discussed by Grebogi, Ott, Pelikan and Yorke |35 



Xj+i = 2a cos 27r6'j tanh Xi (16) 
6,^1 = {io + e,}. (17) 

the reasoning that establishes the existence of strange nonchaotic motion is as follows. For lj 



an irrational number, there are no periodic orbits in this system. The mapping x — 2a tanh x 
is 1-1 and contracting, taking the real line into the interval [—2a, 2a]. The dynamics in 9 
is ergodic in the unit interval and therefore the attractor of the dynamical system must be 
contained in the strip [—2a, 2a] ® [0, 1]. This region contains an invariant subspace, namely 
the line a; = 0. A point Xn with the corresponding = 1/4 will map to {xn+i = 0, On+i = 
uj + 1/4), and hence the subsequent iterates will all remain in this invariant subspace. At the 
same time, the line x = can be made unstable by increasing a: the transverse Lyapunov 
exponent is A = In |a|. Therefore, for |a| > 1 the transverse Lyapunov exponent is positive 
and the line x = is no longer attracting. The total Lyapunov exponent can, however, 
become negative for sufficiently large a so there is an attractor which lies in the region 
X = ±2a, with a dense set of points on the line a; = 0, as well as points off this line. The 
above arguments can be put on firm mathematical footing to show that the attractor is both 
strange and nonchaotic [^ for a > 1. For the Harper map, similar arguments [0,^ can 



be advanced to suggest that for large coupling the dynamics is on SNAs. 



In systems such as the quasiperiodically driven circle map [40| or the logistic map [41 



Xn+i = a[l + ecos27r6'n]a:n(l - x„) (18) 

9n+l = {9n + Uj}, (19) 

SNAs are frequently observed in the neighborhood of the transition to chaos [^,^], so it 



is important to determine the scenarios through which they can be formed, as well as the 
different methods used for their characterization. To date there are several routes known 
for the creation of SNAs; these have been reviewed recently ||3^ . SNAs can be characterized 
by calculating the Lyapunov exponents, fractal dimensions, as well as correlation functions 
and related quantities [02 . 



For map of Harper type, namely of the form Eq. (|1^-|T^) or Eq. (|15]), the nontrivial 
Lyapunov exponent can be easily computed as 

N-l 

A = jnn^^lnx?+i (20) 



while the other Lyapunov exponent is trivially zero. Note that since the dynamics is invert- 
ible, there can be no chaos in this system. 



III. THE FIBONACCI CHAIN 

The Fibonacci chain is the simplest example of a quasicrystal in one-dimension, and 
can be thought of as a two-component lattice. Using the two symbols a and 6, with initial 
sequences 5*0 = h,Si = a, the recursive substitution a ^ ab,b ^ a gives 

Si+i = {Si ■ Si^i} (21) 

with {■} representing concatenation. It is clear that the sequence has elements, where 
Fk is a Fibonacci number specified by Fk+i = Fk + F^-i, Fq = 1, Fi = 1. Labeling sites of a 
lattice of length F^ by the element of a sequence Sk, and specifying the site potential Vi = a 
if the symbol is a, Vi = —a if the symbol is b gives the Fibonacci chain of length F^. We 
are interested in the eigenvalue spectrum of this chain in the limit k oo. 

The Fibonacci chain is the particular case of = 7 = (V5 - l)/2, in the more general 
Kohmoto model, Eq. (p!2D. For the Fibonacci chain, and more generally for any irrational 
u in the Kohmoto model, it is known that all states are critical for any a |I6|JT7|]. We 



can therefore use the alternate 'quantization condition', {E : X{E,uj,a) = 0} to obtain the 
eigenvalue spectrum of the Kohmoto model, as shown in Fig. 1 for the value a = 1. The 
Cantor set structure which is also evident in the spectrum for all irrational uj is more clearly 
seen in the dependence of A on ii^ at fixed uj, a, as in Fig. 2. The curve X{E) meets the line 
A = on a Cantor set of points, namely the spectrum, and in the spectral gaps, the curve 
is parabolic in shape. 

The gaps of the spectrum can be studied for any irrational u; we consider here the case 
of the golden mean since this case has been considered in detail earlier. The behavior of gap 
widths as a function of a has been of interest, not just in the Kohmoto model, but in the 



Harper equation and other related models potentials derived from aperiodic sequence 
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as well iT9|J44 



Note that since the Lyapunov exponent is strictly nonpositive in this system, the dy- 
namics is entirely on attractors (strange as well as non-strange). Following Johnson and 



Moser it is possible to define a winding number as a function of as ||45 

fi(E)= hm ^Wix,,...,x^) (22) 

Af— >oo ZI\ 

where W{xi, . . . ,xn) is an indicator function that counts the number of changes of sign in 
the orbit xi, . . . ,xn- Since the variable Xk is the ratio of the wavefunction at neighbouring 
sites, Q{E) essentially counts the number of nodes in the wavefunction per unit length, and 
is therefore the normalized integrated density of states. The winding number is shown as a 
function of E for a = 1 in the Fibonacci chain in Fig. 3; within a spectral gap, the winding 



number which remains constant, can be expressed as p + q^y, with p,q integers [^. This 
provides a labeling of the gaps which are indicated for a few of the larger gaps in Fig. 3. 

For fixed uj, the ordering of the individual eigenstates do not change with the potential 
strength a. However gap widths w vary with a in a complicated manner. The width of the 
largest gap, namely the one marked A in Fig. 2, scales linearly: wa ~ a. Those marked B 
and C asymptote to constant width, wb,wc ~ while all other gaps decrease in width, 
but as power, w ~ the exponent fi depending on the particular gap. The dependence 
of width on a is shown for gaps A — D in Fig. 4. In the case of Harper equation it has 



been shown I^T] that each gap can be indexed by distinct topological integer index, which 
determines the scaling of the widths. Scaling results are also available for the Thue-Morse 
and period-doubling potentials [|l9l. In the present case it has not been possible to deduce 
any simple number-theoretic dependence of the gap exponents as in the Harper case, though 
clearly, here too the gaps never close since the spectrum is always singular continuous. 

The depths characterising the gaps, namely the minimum value that the LE takes inside 
a gap also vary with a in a similar fashion, increasing for A, asymptoting to constant for B 
and C, and decreasing as a power for all other gaps. 
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IV. SUMMARY AND DISCUSSION 



The equivalence between critical SNAs and critically localized states provides a new 
means of determining the spectrum of discrete quasiperiodic systems. We have, in previous 



studies |]3T| , |43[| , extensively applied this technique to the study of the Harper equation, 
and here we make application to the case of the Kohmoto model, specifically the case of the 
Fibonacci chain, where the potential takes values ±a in a quasiperiodic manner as described 
in Section fT\\ . 

Similar dichotomous potentials which have been studied in the context of critical localiza- 
tion are the Thue-Morse, period-doubling and Rudin-Shapiro sequences ||19[. The methods 
employed here, namely the derivation of an iterative mapping equivalent to the discrete 
Schrodinger equation and the identification of critically localized states with critical SNAs 
can be applied to this entire class of problems. 

The present method makes it possible to empirically study the eigenvalue spectrum of 
such systems in considerable detail. As in related quasiperiodic problems the scaling 



of the gap widths with potential strength a depends on the particular gap being considered. 
For the Fibonacci chain the scaling is w ~ a'^, with the exponent /i = 1 for the largest 
gap, /i = for the next two, and /i < — 1 for the remaining gaps. The organization of the 
gaps, which can be uniquely labeled by two integers, remains incompletely understood, but 



it is likely that, as in the Harper problem [31], the exponent governing the power-law decay 



of gap widths will be related to the gap labels. These aspects, and the application of the 
present technique to more general aperiodic sequences is the subject of ongoing work . 
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Figure Captions 
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FIGURES 
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FIG. 1. Phase diagram for the Kohmoto model for a = 1. 
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FIG. 2. Lyapunov exponent versus energy at a = 1 for the Fibonacci chain, namely 
u = (\/5 — l)/2. The largest visible gaps are labelled A,B,C and D respectively. The dynamics 
for A = corresponds to SNAs. 
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E 

FIG. 3. Integrated density of states as a function of energy E fox lv = j = {^/5-l)/2 at a = 1. 
On the plateaux corresponding to the gaps A,B,C and D the rotation number can be expressed 
as p + (77, with p, q integers. These are indicated for the largest gaps. 
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FIG. 4. Scaling of the gap widths, w in the Fibonacci chain for the largest gaps as a function 
of a. 
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